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Abstract.—The humpback whale 
(Megaptera novaeangliae) is a cosmo- 
politan species whose stocks were dras- 
tically decreased by commercial whal- 
ing practices prior to 1967. The North 
Pacific population was estimated to be 
between 15,000 and 20,000 animals 
before the practice of whaling. At the 
time of the commencement of its inter- 
national protection in 1967, this popu- 
lation may have been reduced to fewer 
than 1000 individuals, The Pacific coast 
of Mexico and the Revillagigedo Archi- 
pelago constitute one of the main breed- 
ing and calving areas for North Pacific 
humpback whales. The objective of this 
paper is to present an estimation of 
abundance of humpback whales in this 
region based on photographic identifi- 
cation of individual animals. Estimates 
of population size were obtained by us- 
ing mark and recapture models for both 
closed and open populations, with each 
year representing a capture occasion. 
A total of 1184 humpback whales were 
identified in Mexican waters between 
1986 and 1993. The best estimates of 
population size for the Mexican stocks 
were those provided by the modified 
Jolly-Seber method: 1813 (95% CI: 918- 
2505) for the coastal stock in 1992, and 
914 (95% CI: 590-1193) for the Revil- 
lagigedo stock in 1991. 
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Humpback whales, Megaptera no- 
vaeangliae, make seasonal migra- 
tions between low-latitude winter- 
ing areas used for mating and calv- 
ing and high-latitude feeding areas. 
The general distribution of feeding 
areas in the North Pacific covers 
coastal waters in the western North 
Pacific from northern Japan 
throughout the Bering Sea and in 
the eastern North Pacific as far 


south as southern California. Dur- 
ing the winter breeding season, 
these whales congregate in three 
geographically isolated tropical ar- 
eas: the Ryukyuan, Bonin, and 
Mariana Islands south of Japan; the 
islands of the Hawaiian Archi- 
pelago; and the Pacific coast of 
Mexico and the Revillagigedo Archi- 
pelago (Rice, 1974; Johnson and 
Wolman, 1984). 


1018 


Fishery Bulletin 97(4), 1999 





Humpback whales have been listed as endangered 
since severe reduction of all stocks worldwide by com- 
mercial exploitation (Rice, 1974; Gambel, 1976). The 
number of these whales were estimated to be between 
15,000 and 20,000 animals before whaling depleted 
them during the first half of the 20th century; at the 
time of its international protection in 1967 this popu- 
lation may have been reduced to fewer than 1000 
individuals (Rice, 1974, 1978). However, the reliabil- 
ity of these figures is unknown. 

The extent of the recovery in the North Pacific 
population over the last 25 years is debatable. Esti- 
mates of abundance with mark and recapture tech- 
niques based upon photo-identification data have 
been made for different areas of this ocean, but there 
has been much debate regarding the reliability of 
such estimates (Darling et al., 1983; Baker et al., 1986; 
Darling and Morowitz, 1986; Baker and Herman, 1987; 
Alvarez et al., 1990; Calambokidis et al., 1990; Cerchio, 
1998; Calambokidis et al.!). There is currently insuffi- 
cient evidence to assess whether a significant change 
in abundance has occurred since whaling ceased. 

Three main wintering aggregations of humpback 
whales are recognized off the Pacific coast of Mexico: 
the Baja California Peninsula; the mainland coast 
of Mexico (including, Isabel Island, Tres Marias Is- 
lands, and the mainland coast); and the Revillagigedo 
Archipelago (including Socorro, San Benedicto, Roca 
Partida, and Clarion Islands) (Rice, 1979; Urban y 
Aguayo, 1987). 

The comparison of photo-identified humpback 
whales in these wintering grounds showed that there 
was a greater affinity between whales off Baja Cali- 
fornia and those off the mainland coast of Mexico 
than those off either Baja California or the main- 
land coast and those off the Revillagigedo Archipelago 
(Urban et al., 1989; Ladrón de Guevara et al., 1993; 
Jaramillo, 1995) (Table 1). From these results, two 
different population units of humpback whales dur- 
ing winter in Mexican waters were previously pro- 
posed: the coastal Stock (including Baja California 
and mainland coast of Mexico), and the Revillagigedo 
stock (Alvarez et al., 1990; Urban et al.?) (Fig. 1). 


Calambokidis, J., G. H. Steiger, J. Straley, T.J Quinn H, L.M 
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structure of humpback whales in the North Pacifie basin. 
Final Rep. to Southwest Fisheries’Science Center, Natl. Mar 
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Table 1 


Matches (above the diagonal) and interchange index’ (be- 
low the diagonal) among the three main wintering aggrega- 
tions in the Mexican Pacific. Sample size in parentheses. 





Baja Mainland Revilla- 
California coastof gigedo 
Peninsula Mexico Archipelago 


(471) (383) (450) 


Baja California Peninsula — 64 
Mainland coast of Mexico 0.38 — 
Revillagigedo Archipelago 0.12 0.18 
1 Index of interchange. 
This index quantify the degree of interchange among two samples 


(among regions) that accounted for sample size: 


Index of interchange = (mi (ny * ng) x 1000, 


where n, = whales identified (captured) in sample 1, 
na = whales identified in sample 2; and 
m,= captured whales from sample 1 recaptured in sample 2. 


(see Baker et al. 1985; Cerchio et al., 1998). 





This division of winter aggregations is also supported 
by mitochondrial DNA lineage analyses (Medrano- 
Gonzalez et al., 1994, 1995a, 1995b); however, be- 
cause the waters off California, Oregon and Wash- 
ington are the primary feeding destination of the 
whales observed in Baja California and the main- 
land coast, the migratory destination of the whales 
seen around the Revillagigedo Islands is still un- 
known (Urban et al. 1987; Urban et al.*). 

Previous estimates of humpback whale abundance 
in Mexican waters are the following: 500-600 for 
Socorro Is. (Campos, 1987); 200-400 for Isabel Is. 
(Alvarez, 1987); 600-700 for the mainland coast 
(Alvarez et al., 1990); and 1200-1700 for all the Mexi- 
can Pacific (Urban et al., 1989). However, it 1s recog- 
nized that these estimates were based on limited data 
and on assumptions that generally were not tested. 

In our study we present an analysis of photographic 
data obtained during the winter breeding and calv- 
ing seasons from 1983 to 1993 in the different as- 
sembly areas of the Mexican Pacific. We use this 
analysis to calculate reliable independent estimates 
of abundance of humpback whales for the coastal and 
the Revillagigedo stocks. 


‘Urban R.J., A. Jaramillo L., A. Aguayo L., Paloma Ladron de 
Guevara P, M. Salinas Z., C. Alvarez F., L. Medrano G., J. 
Jacobsen, K. C. Balcomb, D E. Claridge, J. Calambokidis, G. 
H. Steiger, J. M. Straley, O von Ziegesar, S. Mizroch, M. 
Dahlheim, J. M. Waite, J. D. Darling, and C. S. Baker. 
19xx. Migratory destination of the Mexican Pacific humpback 
whales. Universidad Autonoma de Baja California Sur, Ap. 
Post. 19-B, La Paz, Baja California Sur, Mexico. 
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Figure 1 
Study area showing, by shading, the distribution of the Revillagigedo and coastal stocks. 


Materials and methods 


Individual identification 


Humpback whales were individually identified from 
photographs taken by different institutions from 
approximately December to March, along the main- 
land coast of Mexico (1983-92), from mid-January 
to March in Baja California (1987-93), and at the 
Revillagigedo Archipelago (1986-92) (Table 2). 

Although photo-identification data have been 
available for the coastal stock since 1983 and for the 
Revillagigedo stock since 1985, efforts before 1986 
were limited in time and space, and also lacked con- 
tinuity. Beginning in 1986, surveys were conducted 
off Socorro Island for the Revillagigedo stock, and 
near San José del Cabo or the mainland coast (or 
both) for the coastal stock. No data were obtained 
for either Revillagigedo or the mainland in 1993. 

The whales were identified individually through 
photographs of the black and white pigmentation pat- 
terns on the ventral surface of their flukes (see Katona 
and Whitehead, 1981). Photographs were obtained with 
35-mm cameras equipped with 200-300 mm lenses. 
Although the film used varied by region and season, 
most photographs were taken with black and white 
Kodak Tmax 400 ISO pushed to 1600 ISO, ensuring 
shutter speeds as high as 1/1000 of a second. 


On the basis of focus, angle, and light conditions, 
all fluke photographs were judged to be either good, 
fair, or of poor quality; only photos in the first two 
categories were included in this study. Within these 
quality levels, the whales showed at least 50% of each 
fluke at a sufficiently vertical angle to allow the shape 
of the trailing edge to be distinguished. Calves were 
excluded from the analysis owing to their tendency 
to show changes in pigmentation patterns during the 
first year of life (Carlson et al., 1990). Selected pho- 
tographs were compared visually by at least three 
persons with experience in matching humpback 
whales flukes photographs from both Universidad 
Nacional Autonoma de Mexico and Universidad 
Autonoma de Baja California Sur. 


Abundance estimation 


Abundance estimates were obtained by using an 
eight-year period for the coastal stock, from 1986 to 
1993, and a seven-year period for the Revillagigedo 
stock, from 1986 to 1992. 

Estimates of population size were obtained by us- 
ing mark and recapture models for both closed and 
open populations, with each year representing a cap- 
ture occasion. The time span of the study was seven 
years; consequently, although calves were not in- 
cluded in the analysis, it is inevitable that additions 
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Table 2 
Study periods and data sources for each of the regions. UNAM = Universidad Nacional Autonoma de Mexico; UABCS = Universidad 


Autonoma de Baja California Sur; CWR = Center for Whale Research; CS =Cousteau Society. 





Region Season 


Source 


Study period 





1986 
1987 
1988 
1989 
1990 
199 
1992 


Baja California Peninsula 1987 
1988 
1989 
1990 
1991 
1992 
1993 


1986 
1987 
1988 
1989 
1990 
1991 
1992 


Revillagigedo Archipelago 


to the population were occurring in the form of ani- 
mals born in previous years that became identifiable 
as they grew. Also, some individuals died. Because 
we would therefore expect to be studying an open 
population, we chose the Jolly-Seber mark and re- 
capture model. One of the critical assumptions of this 
method is the highly unlikely condition that all ani- 
mals have the same capture probabilities at the mo- 
ment of the sample (Seber, 1982). If heterogeneity in 
capture probabilities is present, the open population 
models will underestimate population size to a 
greater extent than closed population models (Car- 
others, 1973; Pollock et al., 1990). Given that there 
was no way to account for such problems with open 
population models (Hammond, 1986), we decided to 
use closed population models to test for variations 
in capture probabilities. 

Closed population estimates were obtained and 
assumptions were tested by using the software pro- 
gram CAPTURE developed by Otis et al. (1978), 
which included an algorithm to select the appropri- 
ate model after the hypothesis testing procedure. The 
conceptual basis for this selection procedure is a 
tradeoff between precision and bias. If a simple 
model, such as the null model M, in Otis et al. (1978), 
is used to estimate parameters from data that vio- 


20 Dee 85-4 Mar 86 
25 Dec 86-20 Mar 87 
22 Jan 88-8 Feb 88 
24 Jan 89-5 Mar 89 
25 Nov 89-21 Mar 90 
28 Nov 90-26 Feb 91 
26 Dec 91-6 Mar 92 


UNAM 
UNAM 
UNAM 
UNAM 
CWR, UABCS, UNAM 
UNAM 
UNAM 


10 Feb 87-5 Mar 87 
6 Feb 88-8 Mar 88 
23 Jan 89-25 Mar 89 
20 Feb 90-31 Mar 90 
24 Jan 91-23 Mar 91 
20 Jan 92-2 Apr 92 
15 Jan 93-5 Apr 93 


CWR, UABCS 
CWR, UABCS 
CWR, UABCS 
CWR, UABCS 
UABCS 
UABCS 
UABCS 


16 Jan 86-20 Feb 86 
20 Jan 87-5 Mar 87 

1 Feb 88-9 Mar 88 
16 Jan 89-7 Mar 89 
31 Jan 90-20 Mar 90 
15 Jan 91-25 Apr 91 
25 Apr 92-7 May 92 


UNAM 
UNAM 
UNAM 
UNAM 
UNAM 
UNAM 
UNAM, UABCS, CS 





lates in any way the assumption of equal capture 
probability, then significant biases are introduced in 
parameter estimates and sampling variances will be 
artificially small. On the other hand, if a more com- 
plex model is used, such as M,, of Otis et al.. (1978), 
that allows capture probabilities to vary with time 
and among individuals, biases may be reduced but 
the sampling variance will be greater than it should 
be. The selection procedure takes into account the 
individual goodness-of-fit tests performed for specific 
models on the data and the confrontation of related 
models (i.e. where one model is a particular case ofa 
general one). The significance levels for all these tests 
are combined in a standard discriminant analysis 
and the resulting statistic is standardized so that its 
value ranges from 0 to 1, 1 being the score that indi- 
cates the appropriate model (Otis et al., 1978). 

The basic or null model M, can be applied to the 
general case of multiple recaptures in a closed popu- 
lation, where all animals have equal capture prob- 
abilities, and where this probability remains constant 
in time. The model that allows the capture probabil- 
ity to vary in time, simultaneously permitting indi- 
viduals to have unequal capture probahility (Chao 
et al., 1991: model M,, in Otis et al., 1978) was se- 
lected for estimation of both stock sizes. In addition, 
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abundance was also estimated for the 
coastal stock by using the model that only 
allows relaxation of the requirement for 
constant capture probability in time 
(Darroch, 1958: mode! M, in Otis et al., 
1978). Model Mp provides the estimator 
of population size for a situation where 
capture probabilities vary in time and 
among individuals. Mode! M, provides the 
estimate of population size when capture 
probabilities are the same among individu- 
als but vary in time (Otis et al., 1978). 
Outputs from CAPTURE were used only 
as a “screening technique,” as suggested by 
Menkins and Anderson (1988) to investigate 
departures from the assumption of equal catchability. 
The test for population closure within CAPTURE was 
ignored because its power is low (Otis et al., 1978). 
Open population estimation was done through the 
software program RECAP? This program provides 
estimates of parameters under the basic Jolly-Seber 
model where all individuals have equal capture prob- 
abilities and survivorship, but these were allowed to 
vary between sampling occasions. Also, it incorpo- 
rates a modification of the Jolly-Seber model (modi- 
fied J-S model) that constrains estimates to feasible 
values, stabilizing them and providing more reliable 
confidence intervals (Buckland, 1980). Open popu- 
lation estimates were also obtained by using the soft- 
ware program JOLLY described by Pollock et al. 
(1990). Of particular interest is the goodness-of-fit 
tests performed for this program to investigate how 
well our data are explained by the Jolly-Seber model 
or any of the variants included in the program. The 
results obtained under such variants were also very 
useful because they allowed for the estimation of 
more precise parameters when either capture prob- 
abilities or survivorship (or both) were kept constant. 
Even in the case when the more general Jolly-Seber 
mode} was chosen (model A), observation of survi- 
vorship estimates under model] 2 (temporary trap 
response model) were useful to look for the amount 
of transit within the areas occupied by each stock. 


Results and discussion 


Photo identification 


A total of 1184 humpback whales were identified in 
Mexican waters between 1986 and 1993. Of the to- 


4 Contact Stephen Buckland, Mathematical Institute, North 
Haugh, University of St Andrews, St. Andrews, Fife KY16955, 
U.K. 


Region 


Coast 
Animals caught 4 83 
Newly caught Bf 3 5 91 








Table 3 


Number of hnmpback whales identified in Mexican waters from 1986 
to 1993. 


1986 1987 1988 1989 1990 1991 


1992 1993 


109 


Revillagigedo 
Animals caught 
Newly caught 


Table 4 


Estimates of population size (N) with closed population 
models (see “Materials and methods” section for descrip- 
tions of the models). 


Capture prob. P 


Model 95% CI 
Coastal stock 

2039 
1828 
2188 


1728-2194 
1660-2076 
1861-2612 


0.05 

MY 0.06 average 
M 

Revillagigedo stock 

1078 

874 

1308 


0.05 average 


807-1169 
762-1020 
1043-1684 


0.0687 
0.084 average 


0.057 average 


tal, 733 individuals were observed in the area occu- 
pied by the coastal stock, 412 in Revillagigedo, and 39 
in both areas (a summary of the number of individuals 
identified in each area is provided in Table 3). 


Abundance estimation 


The closed population estimators showed that het- 
erogeneity and variations in capture probabilities in 
time are present in both stocks. However, the confi- 
dence intervals built for models M, and M,, overlap 
to a greater extent in the coastal estimates than in 
the Revillagigedo estimates, indicating that both 
models M,,, and M, are equally good for the coasta} 
stock (Table 4). These relatively small capture prob- 
abilities estimated with the closed models are consis- 
tent with a positive bias in the population size esti- 
mates because, as will be discussed, the populations 
were open through the period of study. Therefore, the 
closed population models were useful in showing that 
the assumption of equal catchability is not met. 
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Table 5 


Estimates of population size (N) with the program JOLLY, 
mode] A (general) (see “Materials and methods” section for 





Table 6 


Estimates of population size (N) with the program RECAP, 
modified Jolly-Seber model (sce “Materials and methods” 



























descriptions of the models). 


7° Goodness-of-fit model 2 


Year N 95% CI Capture prob. P 
Coastal stock! 
1987 Zal 119-1424 0.08 
1988 484 215-753 0.14 
1989 mod 382-1220 0.10 
1990 1247 682-1811 0.17 
1991 1660 610-2712 0.06 
1992 1702 488-2917 0.10 
Average 1109 793-1425 0.11 
Revillagigedo stock? 
1987 134 13-256 025 
1988 367 125-609 0.18 
1989 593 288-897 0.17 
1990 456 179-733 0.10 
1991 1569 -140-3277 ONS 
Average 623 268-980 0.16 
l 4° Goodness-of-fit model A Say) P = 0.8754 





section for descriptions of the models). 





Year 


Coastal stock 
1987 
1988 
1989 
1990 
1991 
1992 
Average 


Revillagigedo stock 
1987 
1988 
1989 
1990 
1991 
Average 


425-1142 
327-931 
519-1229 
824-1613 
1168-2686 
918-2505 
889-1406 


49-484 
215-691 
388-825 
359-891 
590-1193 
373-677 


95% CI Capture prob. P 


0.09 
0.13 
0.10 
0.19 
0.06 
0.10 





* 4° Goodness-of-fit model A 


P = 0,0667 


5.968 
TIZ P= Oesge 
13.23 
0.448 





X? Goodness-of-fit model 2 P= 0799 






The Jolly-Seber analysis indicated that the open 
population model explained in a better way the data 
for both Revillagigedo and coastal stocks (Table 5). 
In the Revillagigedo stock the goodness of fit y? test 
on model A (general model) did not reject the null 
hypothesis (P=0.0667), but the result showed that 
data are better explained (P=0.799) by model 2 (tem- 
porary trap response model), assuming that survival 
rates of newly captured animals differ from those of 
previously captured animals. In other words, this test 
indicates that some animals may have moved into 
the capture site and once captured, they could have 
left and even re-entered the sampling area. This in- 
ference an be supported by the fact that data for this 
stock were collected only in Socorro Island and that 
some animals may have a preference for anther is- 
land, such as Clarion, far away from the sampling 
area (Fig. i). For the coastal stock, the goodness-of- 
fit y? tests indicate that the model A (P=0.8754) and 
model 2 (P=0.8892) are almost equally good in ex- 
plaining the data. Although sampling coverage was 
wider for this stock, allowing a better collection of 
photographs for moving animals, this problem still 
remained, but had a smaller effect on the estimates. 

Having accepted that the Jolly-Seber model ex- 
plained the data well for the two stocks, population- 
size estimates were obtained with the software pro- 


gram RECAP. Results show that the modified esti- 
mator (modified J-S model) provided more stable 
estimates and narrower and more realistic confidence 
intervals (Table 6). 

Because capture probabilities are larger than those 
obtained by the closed population methods, and con- 
sidering that the Jolly-Seber estimates may be nega- 
tively biased given heterogeneity, population size are 
likely overestimated by the closed population methods. 

In hght of the above arguments, we conclude that 
the best estimates of population size for the Mexi- 
ean stocks, using the data at hand, are those pro- 
vided by the modified Jolly-Seber mode) of program 
RECAP. Estimates obtained with closed population 
methods are not recommended when data are pooled 
for the whole extension of the study; however, it may 
be worth obtaining estimates with pairs of years with 
models testing for heterogeneity and capture prob- 
abilities varying with time, although certain infor- 
mation and precision may be sacrificed. 

Estimates for each year show an apparent trend 
towards increase for both stocks. Such a trend was 
considered an artifact of a reduction in the bias 
caused by heterogeneity in capture probabilities af- 
ter sample size increased with time. Consequently, 
we consider that the best estimates for the stock sizes 
were those obtained for the last years and are 1813 
(95% CI: 918-2505) for the coastal stock in 1992, and 
914 (95% C]: 590-1193) for the Revillagigedo stock 
im LOOM 
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